We present the first NIR spectro-interferometry of the LBV 
77 Carinae. The observations were performed with the AMBER 
instrument of the ESO Very Large Telescope Interferometer 
(VLTI) using baselines from 42 to 89 m. The aim of this 
work is to study the wavelength dependence of 77 Car's opti- 
cally thick wind region with a high spatial resolution of 5 mas 
(11 AU) and high spectral resolution. The observations were 
carried out with three 8.2 m Unit Telescopes in the /if -band. 
The raw data are spectrally dispersed interferograms obtained 
with spectral resolutions of 1,500 (MR-K mode) and 12,000 
(HR-K mode). The MR-K observations were performed in 
the wavelength range around both the HeI2.059/iin and the 
Bry2.166/iin emission lines, the HR-K observations only in 
\Q • the Bry line region. The spectrally dispersed AMBER interfer- 
' ograms allow the investigation of the wavelength dependence 

■ of the visibility, differential phase, and closure phase of 77 Car. 

, In the /f-band continuum, a diameter of 4.0+0.2 mas (Gaussian 
Qh- FWHM, fit range 28-89 m baseline length) was measured for 
^t. ' 77 Car's optically thick wind region. If we fit Hillier et al. (2001) 

■ model visibilities to the observed AMBER visibilities, we ob- 
tain 50% encircled-energy diameters of 4.2, 6.5 and 9.6 mas in 

■ the 2.17 pm continuum, the He I, and the Bry emission lines, 
I respectively. In the continuum near the Bry line, an elonga- 

tion along a position angle of 120° + 15° was found, consistent 
[ with previous VLTI/VINCI measurements by van Boekel et al. 
,— ' 1 (2003). We compare the measured visibilities with predictions 
of the radiative transfer model of Hillier et al. (2001), find- 
^ ■ ing good agreement. Furthermore, we discuss the detectability 
I of the hypothetical hot binary companion. For the interpreta- 
• tion of the non-zero differential and closure phases measured 
"pH I within the Bry line, we present a simple geometric model of an 
O ( inclined, latitude-dependent wind zone. Our observations sup- 
' I port theoretical models of anisotropic winds from fast-rotating, 
^ luminous hot stars with enhanced high-velocity mass loss near 
the polar regions. 

£>. ' Key words. Stars: individual: 77 Carinae - Stars: mass-loss, 
emission-line, circumstellar matter, winds, outflows - Infrared: 
stars - Techniques: interferometric, high angular resolution, 
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1. Introduction 

The enigmatic object 77 Car is one of the most lumi- 
nous and most massive (M ~ 100 M ) unstable Luminous 
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Blue Variables suffering from an extremly high mass loss 
(Davidson & Humphreys 1997). Its distance is approximately 
2300+100 pc (Davidson & Humphreys 1997; Davidson et al. 
2001; Smith 2006). r\ Car, which has been subject to a va- 
riety of studies over the last few decades, is surrounded by 
the expanding bipolar Homunculus nebula ejected during the 
Great Eruption in 1843. The inclination of the polar axis of the 
Homunculus nebular with the line-of-sight is ~ 41° with the 
southern pole pointing towards us (Davidson et al. 2001 ; Smith 
2006). The first measurements of structures in the innermost 
sub-arcsecond region of the Homunculus were obtained by 
speckle-interferometric observations (Weigelt & Ebersberger 
1986; Hofmann & Weigelt 1988). These observations revealed 
a central object (component A) plus three compact and sur- 
prisingly bright objects (components B, C, and D) at distances 
ranging from approximately 0.1 to 0.2 ". HST observations of 
the inner 1 region (Weigelt et al. 1995) provided estimates 
of the proper motion of the speckle objects B, C, and D (ve- 
locity ~ 50 km/s; the low velocity suggests that the speckle 
objects are located within the equatorial plane), and follow-up 
HST spectroscopy unveiled their unusual spectrum (Davidson 
et al. 1995). The central object (speckle object A) showed 
broad emission lines, while the narrow emission lines came 
from the speckle objects B, C, and D. Therefore, A is cer- 
tainly the central object while B,C, and D are ejecta. Recent 
observations of r\ Car by Chesneau et al. (2005) using NACO 
and VLTI/MIDI revealed a butterfly-shaped dust environment 
at 3.74 and 4.05 fim and resolved the dusty emission from the 
individual speckle objects with unprecedented angular resolu- 
tion in the NIR. Chesneau et al. also found a large amount of 
corundum dust peaked ~ 1 south-east of the central object. 

Spectroscopic studies of the Homunculus nebula showed 
that the stellar wind of r\ Car is aspherical and latitude- 
dependent, and the polar axes of the wind and the Homunculus 
appear to be aligned (bipolar wind model; Smith et al. 2003). 
Using Balmer line observations obtained with HST/STIS, 
Smith et al. (2003) found a considerable increase of the wind 
velocity from the equator to the pole and that the wind density 
is higher in polar direction (parallel to the Homunculus; PA of 
the axis ~132°; Davidson et al. 2001) than in equatorial direc- 
tion by a factor of ~2. van Boekel et al. (2003) resolved the op- 
tically thick, aspheric wind region with NIR interferometry us- 
ing the VLTI/VINCI instrument. They measured a size of 5 mas 
(50% encircled-energy diameter), an axis ratio of 1.25 ± 0.05, 
and a position angle (PA) of the major axis of 134°±7°, and de- 
rived a mass-loss rate of 1.6x 10~ 3 M /yr~'. The aspheric wind 
can be explained by models for line-driven winds from lumi- 
nous hot stars rotating near their critical speed (e.g., Owocki 
et al. 1996, 1998). The models predict a higher wind speed 
and density along the polar axis than in the equatorial plane. 
In addition, van Boekel et al. showed that the K broad-band 
observations obtained with VINCI are in agreement with the 
predictions from the detailed spectroscopic model by Hillier 
etal. (2001). 

The Hillier et al. (2001, 2006) model was developed to 
explain STIS HST spectra. The luminosity of the primary 
(5 x 10 6 L o ) wa s set by observed IR fluxes (see discussion 
by Davidson & Humphreys 1997) and the known distance of 



2.3 kpc to 77 Car. Any contribution to the IR fluxes by a bi- 
nary companion was neglected. Modeling of the spectra was 
undertaken using CMFGEN, a non-LTE line blanketed radia- 
tive transfer developed to model stars with extended outflow- 
ing atmospheres (Hillier & Miller 1998). For the modeling of 
77 Carinae, ions of H, He, C, N, O, Na, Mg, Al, Si, S, Ca, Ti, Cr, 
Mn, Fe, Ni, and Co were included. The mass loss was derived 
from the strength of the hydrogen lines and their associated 
electron scattering wings. Due to a degeneracy between the 
mass-loss rate and the He abundance, the H/He helium abun- 
dance ratio could not be derived, but was set at 5: 1 (by number), 
which is similar to that found by Davidson et al. (1986) from 
nebula studies. CNO abundances were found to be consistent 
with those expected for full CNO processing. With the excep- 
tion of Na (which was found to be enhanced by at least a factor 
of 2), the adoption of solar abundances for other metal species 
was found to yield satisfactory fits to the STIS spectra. A more 
recent discussion of the basic model, with particular reference 
to the UV and outer wind, is given by Hillier et al. (2006). 

Because the wind is optically thick, the models are fairly 
insensitive to the radius adopted for the hydrostatic core (i.e., 
the radius at which the velocity becomes subsonic). One ex- 
ception was the He I lines, which decreased in strength as the 
radius increased and, in general, were very sensitive to model 
details. Additional HST STIS observations show that the He I 
lines are strongly variable and blue-shifted throughout most of 
the 5.54-year variability period. These observations cannot be 
explained in the context of a spherical wind model. It now ap- 
pears likely that a large fraction of the He I line emission origi- 
nates in the bow shock and an ionization zone, associated with 
the wind-wind interaction zone in a binary system (Davidson 
et al. 1999; Davidson 2001; Hillier et al. 2006; Nielsen et al. 
2006). Consequently, the hydrostatic radius derived by Hillier 
et al. (2001) is likely to be a factor of 2 to 4 too small. Because 
the wind is so thick, a change in radius will not affect the Bry 
formation region, and it will only have a minor influence on 
the Bry continuum emitting region. If this model is correct, the 
He I emission will be strongly asymmetrical and offset from the 
primary star. 

A variety of observations suggest that the central source 
of 77 Car is a binary. Damineli (1996) first noticed the 5.5- 
year periodicity in the spectroscopic changes of this object (see 
Damineli et al. 1997, 2000; Duncan et al. 1999; Ishibashi et al. 
1999; Davidson et al. 1999, 2000; van Genderen et al. 2003; 
Steiner & Damineli 2004; Whitelock et al. 2004; Corcoran 
2005; Weis et al. 2005). On the other hand, to date, the bi- 
nary nature of the central object in 77 Car and its orbital pa- 
rameters are still a matter of debate (see, e.g., Zanella et al. 
1984; Davidson 1999, 2001; Davidson et al. 1999, 2000, 2005; 
Ishibashi et al. 1999; Smith et al. 2000; Feast et al. 2001; 
Ishibashi 2001; Pittard & Corcoran 2002; Smith et al. 2003; 
Martin et al. 2006). 

The 1997.9 X-ray peak with the subsequent rapid drop 
to a few-month-long minimum was detected by RXTE (see 
Corcoran 2005). Then the first spectra with HST/STIS were ob- 
tained at 1998.0, demonstrating changes in both the central star 
and the aforementioned speckle objects (Davidson et al. 1999; 
Gull et al. 1999). Pittard & Corcoran (2002) demonstrated that 
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the CHANDRA X-ray spectrum can be explained by the wind- 
wind collisions of the primary star (M = 2 x 10 _4 M o /yr _1 at 
500 km/s) and a hot companion (M = 10~ 5 M Q /yr~ l at 3,000 
km/s). Verner et al. (2005) used models calculated with the 
CLOUDY code to demonstrate that during the spectroscopic 
minimum, the excitation of the speckle objects is supported by 
the primary stellar flux, but that the UV flux of a hot companion 
consistent with an 07.5V, 091, or early WN star was probably 
necessary to excite the speckle objects during the broad spec- 
troscopic maximum. 

In this paper we present the first spectro-interferometric 
/if -band observations of 77 Car obtained with the VLTI beam- 
combiner instrument AMBER with medium and high spectral 
resolution and in the projected baseline range from 28 to 89 m. 

The paper is organized as follows: In Sect. 2 we give an 
overview of the AMBER observations of 77 Car and describe the 
data reduction procedure in detail, and in Sect. 3, the analyses 
of the continuum data and the measurements within the Bry 
and He I lines are discussed individually. 



2. AMBER observations and data processing 

AMBER (Petrov et al. 2003, 2006a,b) is the near-infrared 
(/, H, K band) beam-combiner instrument of ESO's Very 
Large Telescope Interferometer, which allows the measure- 
ment of visibilities, differential visibilities, differential phases, 
and closure phases (Petrov et al. 2003; Millour et al. 2006). 
AMBER offers three spectroscopic modes: low (LR mode; 
R=A/AA=75), medium (MR mode; fl=l,500), and high (HR 
mode; /?=12,000) spectral resolutions. The fibers in AMBER 
limit the field-of-view to the diameter of the fibers on the sky 
(~ 60 mas). In AMBER the light is spectrally dispersed using 
a prism or grating. The AMBER detector is a Hawaii array de- 
tector with 512x512 pixels. 

Figure 1 shows two AMBER raw interferograms taken in 
the wavelength range around the Bry line in HR (top) and MR 
(bottom) mode. In the MR data sets, the Doppler-broadened 
Bry line covers ~ 8 spectral channels, whereas in HR mode, 
the line is resolved by ~ 50 spectral channels. 

77 Car was observed with AMBER on 2004 December 26, 
2005 February 25, and 2005 February 26 with the three 8.2 m 
Unit Telescopes UT2, UT3, and UT4. With projected base- 
line lengths up to 89 m, an angular resolution of ~5 mas was 
achieved in the K band. As listed in Table 1, the MR-K obser- 
vations were performed in the wavelength range around both 
the He I 2.059^/m and the Bry 2.166/im emission lines. The 
HR-K observations were only performed in a wavelength range 
around the Bry line. The widths of the wavelength windows of 
the obtained MR-K and HR-K observations are approximately 
0.05 pm and 0.02 pm, respectively. 

For the reduction of the AMBER data, we used version 
2.4 of the amdlib 1 software package. This software uses the 
P2VM ( pixel-to-visibilities matrix) algorithm (Tatulli et al. 
2006) in order to extract complex visibilities for each baseline 
and each spectral channel of an AMBER interferogram. From 
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1 This software package is available from 
http : / /amber . obs . uj f -grenoble . fr 



Fig. 1. Spectrally dispersed VLTI/ AMBER Michelson interferograms 
of 77 Car. The two panels show the spectrally dispersed fringe signal 
(IF) as well as the photometric calibration signals from the three tele- 
scopes (P1-P3) in high (HR, upper panel) and medium spectral reso- 
lution mode (MR, lower panel). In both panels, the bright regions are 
associated with the Doppler-broadened Bry emission line. 



these three complex visibilities, the amplitude and the closure 
phase are derived. While the closure phase is self-calibrating, 
the visibilities have to be corrected for atmospheric and in- 
strumental effects. This is done by dividing the 77 Car visibil- 
ity through the visibility of a calibrator star measured on the 
same night. In order to take the finite size of the calibrator star 
into account, the calibrator visibility is corrected beforehand 
through division by the expected calibrator star visibility (see 
Table 1). In the case of the MR measurement performed on 
2005-02-25, the interferograms recorded on the calibrator con- 
tain only fringes corresponding to the shortest baseline (UT2- 
UT3). Thus, the 77 Car visibility for this night could only be 
calibrated for this shortest baseline. 

Besides the calibrated visibility and the closure phase, the 
spectral dispersion of AMBER also allows us to compute dif- 
ferential observables; namely the differential visibility and the 
differential phase (Petrov et al. 2003, 2006a,b; Millour et al. 
2006). These quantities are particularly valuable, as they pro- 
vide a measure of the spatial extent and spatial offset of the line- 
emitting region with respect to the continuum emission. Since 
the measured complex visibilities are affected by wavelength- 
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Fig. 2. AMBER observables derived from our 77 Car data around the Bry line for three independent measurements (Left: MR, 
2004 December 26, Middle: MR, 2005 February 25, Right: HR, 2005 February 26). The first row shows the continuum-normalized spec- 
tra as extracted from the interferometric channels, followed by the derived calibrated visibilities and the differential visibilities. In the fourth 
and fifth row, the differential phase and the closure phase are presented. In the spectra we mark the wavelength regimes, which we defined as 
continuum for our analysis (shaded regions). The vertical grey line marks the rest-wavelength of Bry (/l vac = 2.1661 /jm; the small correction 
due to the system velocity of -8 km/s (Smith 2004) has been neglected). We show different error bars within each panel: The left error bars 
correspond to the total (including statistical and systematic) error estimated for the continuum wavelength range, and the error bar towards the 
right visualizes the total error for the wavelength range within the line. For the HR 2005-02-26 measurement, data splitting showed that for a 
small wavelength range (hatched areas in the two lower right panels), the differential phase for the longest and middle baseline as well as the 
closure phase become very noisy and are therefore not reliable. Furthermore, the HR differential phase of the longest baseline is noisy at all 
wavelengths. See Sect. 2 for further details. 
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Table 1. Summary of the AMBER ;/ Car observations using the UT2, UT3, and UT4 telescopes. 



Date [UT] 


Time [UT] 


Orbital 


spectral 


line within 


DIP' 


N r 
t; Car 


Calibrator 


» r S 

N c alib . 


Calibrator 




Start 


end 


phase c 


mode 


spectral range 








uniform disk 














[ms] 








diameter [mas] 


2004 Dec. 26 


07:52 


08:16 


0.267 


MR-K 


Bry 


40 


7,500 


HD 93030 


5,000 


0.39° 




08:19 


08:32 


0.267 


MR-K 


He I 2.059 


40 


5,000 


HD 93030 


5,000 


0.39° 


2005 Feb. 25 


04:33 


04:43 


0.298 


MR-K 


Bry 


50 


5,000 


HD 89682 


2,500 


3.08* 




04:55 


05:05 


0.298 


MR-K 


HeI2.059yum 


50 


5,000 


HD 89682 


2,500 


3.08* 


2005 Feb. 26 


08:16 


08:57 


0.299 


HR-K 


Bry 


82 


7,500 


LCar 


2,500 


2.70 c 



Notes - " Uniform disk (UD) diameter estimated using the method described by Dyck et al. (1996). 
b UD diameter taken from the CHARM2 catalog (Richichi et al. 2005). 

c UD diameter of L Car at the time of the AMBER high-resolution observations derived from the limb-darkened diameter ^ld = 2.80 mas at 
the L Car pulsation phase 0=0.0 (Kervella et al. 2004b, 2006) and d VD /d LD = 0.966 (Kervella et al. 2004a). 

Detector integration time per interferogram. 
' The orbital phase was computed assuming a zero point at JD 2 450 800.0 and a period of 2024 days (Corcoran 2005). 
f Number of 77 Car interferograms. 
8 Number of calibrator interferograms. 



dependent atmospheric piston (optical path difference), the pis- 
ton has to be estimated and subtracted. This was done using the 
ammyorick 1 tool (version 0.56). 

Since a large fraction of the interferograms is of low con- 
trast (probably due to vibration; see Malbet et al. 2006), we 
removed a measurement from the data sets if (a) the intensity 
ratio of two of the photometric channel signals is larger than 4 
(a large ratio means that the interferograms are very noisy since 
the signal is very weak in one channel) or (b) it belongs to the 
70 percent of the interferograms with the lowest fringe contrast 
SNR (with the SNR defined as in TatuUi et al. 2006). In or- 
der to optimize the selection for each baseline of the telescope 
triplet, both of these criteria are applied for each telescope pair 
individually. Furthermore, the first 10 frames in each new se- 
quence of recorded interferograms are removed since they are 
degraded by electronic noise. 

Figures 2 and 4 show the spectra as well as the wavelength 
dependence of the visibilities, differential visibilities, differen- 
tial phases, and closure phases derived from the AMBER inter- 
ferograms for the observations around the Bry and He I emis- 
sion lines. The uv coverage of the observations is displayed in 
Fig. 3. 

The 77 Car spectra were corrected for instrumental effects 
and atmospheric absorption through division by the calibrator 
spectrum. For the HR 2005-02-26 measurement, we found that 
the calibrator itself (L Car) shows prominent Bry line absorp- 
tion (see Fig. A. 1). Therefore, we had to remove this stellar line 
by linear interpolation before the spectrum could be used for 
the calibration. The wavelength calibration was done using at- 
mospheric features, as described in more detail in Appendix A. 

In order to test the reliability of our results, we split each 
of the raw data sets into 5 subsets, each containing the same 
number of interferograms. The results obtained with these in- 
dividual subsets allowed us to test that the major features de- 
tected in the visibility, differential visibility, differential phase, 
and closure phase are stable, even without any frame selection 
applied. As an exception, we found that for a small wavelength 
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Fig. 3. uv coverage of the AMBER rj Car observations. The data ob- 
tained with medium spectral resolution in Dec. 2004 and Feb. 2005 
are indicated by red squares and green bullets, respectively, while the 
high-resolution measurements are shown as blue triangles. Filled sym- 
bols denote observations around the Bry line, and open symbols de- 
note observations around the He I line. 



range of the HR 2005-02-26 data set (hatched areas in the two 
lower right panels of Figure 2), the differential phase corre- 
sponding to the middle and longest baselines and the closure 
phase vary strongly within the subsets and are, therefore, un- 
reliable. This is likely due to the very low visibility value on 
these two baselines, resulting in a low fringe SNR within this 
wavelength range. Furthermore, with this method we found that 
the differential visibility, differential phase, and closure phase 
extracted from the MR 2005-02-25 He I data set are very noisy 
and not reliable. Therefore, these differential quantities and clo- 
sure phases were dropped from our further analysis. 

The subsets were also used to compute statistical errors. 
We estimated the variance for each spectral channel and de- 
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rived formal statistical errors for both the continuum and line 
wavelength ranges. In each panel of Figures 2 and 4, we show 
two types of error bars corresponding to these regions, which 
not only take these statistical errors but also a systematic error 
(e.g. resulting from an imperfect calibration) into account. 



3. Observational results and interpretation 

3.1. Comparison of the observed wavelength 
dependence of the visibility with the NLTE 
radiative transfer model of Hillier et al. (2001) 

For the analysis presented in this chapter, we used the AMBER 
data sets from 2004 Dec. 26 and 2005 Feb. 25 and 26, presented 
in Figs. 2 and 4, and compared the AMBER visibilities and 
spectra with the NLTE radiative transfer model of Hillier et al. 
(2001). To directly compare the AMBER measurements with 
this model, we derived monochromatic model visibilities for all 
wavelengths between 2.03 and 2.18/mi (with AA = 10~ 4 /iin) 
from the model intensity profiles, assuming a distance of 2.3 
kpc for j] Car. The comparison is visualized in Fig. 5 for the 
individual AMBER HR and MR measurements. The first row 
displays the AMBER and model spectra, while all other pan- 
els show the AMBER and model visibilities for the different 
projected baselines. We note that for the comparison shown 
in Fig. 5, we used the original model of Hillier et al. (2001) 
without any additional size scaling or addition of a background 
component. 

As the figure reveals, the NLTE model of Hillier et al. 
(2001) can approximately reproduce the AMBER continuum 
observations for all wavelengths (i.e. 2.03-2.18/mi) and all 
baselines. Moreover, the wavelength dependence of the model 
visibilities inside the Bry line is also similar to the AMBER 
data. There is a slight tendency for the model visibilities in the 
Bry line to be systematically lower, which can be attributed to 
the overestimated model flux in the line. On the other hand, 
there is an obvious difference in the wavelength dependence 
of the visibility across the He I line between the observations 
and the model predictions. This difference probably indicates 
that the primary wind model does not completely describe the 
physical origin and, hence, the spatial scale of the He I line- 
forming region. The discrepancy is possibly caused by addi- 
tional He I emission from the wind-wind interaction zone be- 
tween the binary components and by the primary's ionized 
wind zone caused by the secondary's UV light illuminating the 
primary's wind (e.g., Davidson et al. 1999; Davidson 2001; 
Pittard & Corcoran 2002; Steiner & Damineli 2004; Hillier 
et al. 2006; Nielsen et al. 2006; Martin et al. 2006), as discussed 
in Sects. 1, 3.4.2, and 3.7.2 in more detail. 

Figure 6 shows the AMBER and model visibilities as a 
function of spatial frequency and the corresponding model 
center-to-limb intensity variations (CLVs) for seven selected 
wavelengths (2 continuum wavelengths; center, blue-shifted, 
and red-shifted wings of Bry emission; center of both He I 
emission and absorption). As Fig. 6 reveals, at several wave- 
lengths we find a very good agreement between the visibilities 
measured with AMBER and the visibilities predicted by the 
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Fig. 5. Comparison of the AMBER spectra and visibilities with the NLTE model predictions of Hillier et al. (2001). The figure displays the 
spectra (upper row) and visibilities (lower three rows, see labels for projected baselines) of the four AMBER measurements (green lines) and 
the corresponding data of the Hillier et al. NLTE model (red lines). The errors of the AMBER continuum and line visibility measurements 
are indicated by the two vertical error bars (see Figs. 2 and 4; the left bar is the continuum error bar), and the uncertainty of the AMBER 
wavelength calibration is indicated by the horizontal error bar. As the figure shows, we find good agreement between the AMBER data and the 
model predictions for the continuum visibilities as well as the shape and depth of the visibility inside the Bry line. In the case of the He I line, 
the wavelength dependence of the model visibility inside the line differs considerably from the AMBER measurements, indicating a different 
physical process involved in the line formation (see Sect. 3.1). The Hel wavelength shift can be attributed to a combination of both Doppler 
shift (e.g. Nielsen et al. 2006; Hillier et al. 2006) and uncertainties in the wavelength calibration of the AMBER data. Note that no additional 
scaling has been applied to the Hillier et al. model. The model spectra and visibilities have a spectral resolution (A/AA ~ 20, 000) comparable 
to the HR measurements. 



model of Hillier et al. (2001). This is especially true for the 
continuum data (upper two panels). 

^From the model CLVs, FWHM model continuum diame- 
ters of 2.24 mas and 2.33 mas can be derived for A - 2.040 and 
2. 174 /mi, respectively. If we allow for a moderate rescaling of 
the size of the model, we find that the best x 2 fit at both contin- 
uum wavelengths can be obtained with scaling factors of 1.015 
and 1 .00, respectively. This means that the model size has to be 
increased by only 1.5% at A — 2.040pm and that the best fit 
at 2. 174 /mi is indeed obtained with the original Hillier model 
with a scaling factor of 1 .0. Thus, taking the slight rescaling for 
the best^f 2 fit into account, we can conclude that, based on the 
NLTE model from Hillier et al. (2001, 2006) and the AMBER 
measurements, the apparent FWHM diameters of r\ Car in the 



/if -band continuum at A = 2.040 fim and 2. 174 //m are 2.27 mas 
and 2.33 mas, respectively (see Table 2), corresponding to a 
physical size of approximately 5 AU. 

Since the deviations between the model and the measure- 
ments are larger in the case of the Bry and He I line data (lower 
5 panels in Fig. 6), the scaling factors corresponding to the best 
X 1 fit in the lines show stronger deviations from unity. For the 
Bry emission line, we find scaling factors of 0.74, 0.76, and 
0.78 for A = 2.1661, 2.1635, and 2. 1669 /mi, corresponding to 
FWHM diameters of 1.83, 9.52 and 2.02 mas (see Table 2). 

For the Hel emission line, rescaled models with scaling 
factors of 1 . 24 and 1.11 provide the best x 1 fit for the peaks of 
the emission and absorption within the He I line (A = 2.058 and 
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Fig. 6. Left: Comparison of the AMBER 
visibilities (filled green squares; baseline 
range 28-89 m) as a function of spatial fre- 
quency with the NLTE model predictions of 
Hillier et al. (2001) (solid red lines) for two 
continuum wavelengths (upper two panels; 
see labels for the exact wavelengths), the 
central wavelength of the Bry emission line 
(third), wavelengths in the blue and red 
wing of the Bry emission line (fourth and 
fifth row; at these wavelengths, the AMBER 
data show the strongest differential and clo- 
sure phase signals), and the central wave- 
lengths of the Hel emission and P Cygni 
absorption (lower two panels). The blue tri- 
angles are the background-corrected VINCI 
A'-band measurements from van Boekel 
et al. (2003). Right: Center-to-limb varia- 
tion (CLV; i.e. intensity as a function of an- 
gular radius) of the monochromatic Hillier 
et al. (2001) NLTE models for the wave- 
lengths indicated by the labels. As the fig- 
ures show, the typical FWHM diameter of 
the models is of the order of ~2^1 mas for 
an assumed distance of 2.3 kpc for rj Car. 
See text for details. 



2.056 jum), resulting in FWHM diameters of 4.24 and 4. 19 mas, 
respectively. 

In addition to the inner CLV core, at several wavelengths, 
the CLVs show a very extended wing corresponding to the ex- 
tended Bry and He I line emission regions. Since the intensity 
in the wing is much lower than 50% of the peak intensity, the 
FWHM diameter is not very sensitive to this part of the CLV. 
In other words, in the case of CLVs with multiple or very ex- 
tended components, a FWHM diameter can be quite mislead- 
ing. In such a case, it seems to be more appropriate to use, for 



instance, the diameter measured at 10% of the peak intensity 
(a?io%) or the 50% encircled- energy diameter (^5o%eed)- For 
example, at A — 2.1661 /mi we obtain d\o% = 9.39 mas and 
dso% eed = 9.58 mas, while for the continuum at 2. 174 /mi we 
find dio% = 5.15 mas and d$o% eed = 4.23 mas. Thus, based 
on <f 50% eed, t] Car appears ~ 2.2 times larger at A — 2.1661 /mi 
compared to the continuum at A = 2.17 '4 /mi. The best-fit model 
diameters at the other wavelengths are listed in Table 2. The er- 
rors of the diameter measurements are +4% for the two contin- 
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Table 2. Diameters obtained by fitting Hillier et al. (2001) model vis- 
ibilities to the measured AMBER visibilities. Errors are ±4% for the 
diameters in the continuum and ±10% in the lines (see text). 



Spectral region 


Wavelength 


^FWHM 




^50%EED 




[fim] 


[mas] 


[lllcis] 


[mas] 


continuum 


2.0400 


2.27 


4.85 


3.74 


continuum 


2.1740 


2.33 


5.15 


4.23 


Bry (center) 


2.1661 


1.83 


9.39 


9.58 


Bry (blue wing) 


2.1635 


9.52 


16.46 


9.60 


Bry (red wing) 


2.1669 


2.02 


9.61 


9.78 


He I (absorption) 


2.0560 


4.24 


8.22 


5.36 


He I (emission) 


2.0580 


4.19 


4.30 


6.53 


^fwhm = FWHM diameter; rf 10 % = 


diameter measured at 10% p> 



intensity; c/50% eed = 50% encircled-energy diameter. 

uum diameters and +10% for the line diameters, derived from 
the visibility errors and the uncertainty of the fitting procedure. 

3.2. Continuum visibilities 

3.2.1 . Comparison of the continuum visibilities with the 
Hillier et al. (2001) model predictions 

The comparison of the AMBER continuum visibilities with the 
NLTE model from Hillier et al. (2001, 2006) is shown in the 
two upper left panels of Fig. 6 for the continuum near the He I 
2.059 /im and Bry 2.166/im emission lines (the exact wave- 
lengths are described in Fig. 6). Taking a slight rescaling into 
account, we concluded in the previous section that, based on 
the NLTE model from Hillier et al. (2001) and the AMBER 
measurements, the apparent 50% encircled-energy diameters 
^50% eed of r\ Car in the /if -band continuum at A = 2.040 /an and 
2.174//m are 3.74 mas and 4.23 mas, respectively (see Table 2). 
These diameters are in good agreement with the 50% encircled- 
energy K-band diameter of 5 mas reported by van Boekel et al. 
(2003). 

For comparison, we also fitted the AMBER visibilities with 
simple analytical models such as Gaussian profiles, as de- 
scribed in more detail in Sect. B in the Appendix. From a 
Gaussian fit of the AMBER visibilities, we obtain a FWHM 
diameter of c/Gauss ~ 4-0 ± 0.2 mas in the /f-band continuum. 
As outlined in Sect. B, the diameter value strongly depends 
on the range of projected baselines used for the fit, since a 
Gaussian without an additional background component is not 
a good representation of the visibility measured with AMBER. 
As discussed in Appendix B, using a Gaussian fit with a fully 
resolved background component as a free parameter results in 
a best fit with a 30% background flux contribution (see also 
Petrov et al. 2006a). 

3.2.2. Comparison of the VINCI and AMBER 
continuum visibilities 

In Fig. 6 (left, second row) displaying the averaged Bry contin- 
uum data, the visibilities of 77 Car obtained with VLTI/VINCI 
are shown in addition to the AMBER data. These VINCI mea- 



surements were carried out in 2002 and 2003 using the 35 cm 
test siderostats at the VLTI with baselines ranging from 8 to 
62 m (for details, see van Boekel et al. 2003). Like AMBER, 
VINCI is a single-mode fiber instrument. Therefore, its field- 
of-view is approximately equal to the Airy disk of the tele- 
scope aperture on the sky, which is ~ 1.4" in the case of the 
siderostats. From the VINCI measurements and using only the 
24 m baseline data, van Boekel et al. (2003) derived a FWHM 
Gaussian diameter of c/Gauss ~ 7 mas for the wind region of 
77 Car. At first glance, this diameter measurement seems to con- 
tradict the c/Gauss ~ 4.0 mas FWHM diameter derived from the 
AMBER data. This is not the case, however, since the diam- 
eter fit is very sensitive to the baseline (or spatial frequency) 
fit range, because a Gaussian is not a good representation of 
the visibility curve at all, as can be seen in Fig. B.l. If only 
the VINCI data points are fitted, which have spatial frequen- 
cies < 60 cycles/arcsec (corresponding to projected baselines 
< 28 m), c/Gauss ~ 7 mas provides the best fit. On the other 
hand, if the data point at 136 cycles/arcsec (corresponding to a 
projected baseline of ~ 62 m) is included in the fit, we obtain 
tfcauss ~ 4.3 mas (see also the discussion Sect. B). Thus, when 
using comparable baseline ranges for the Gaussian fits, there 
is good agreement between the AMBER and VINCI measure- 
ments. 

To account for the background contamination of the VINCI 
data caused by nebulosity within VINCI's large 1.4" field-of- 
view (in which, for instance, all speckle objects B, C, and D are 
located), van Boekel et al. introduced a background component 
(derived from NACO data) providing 55% of the total flux. 
Adding this background component to the model of Hillier 
et al. (2001), they found a good match between the model and 
the observations. Since our AMBER observations were carried 
out with the 8.2 m Unit Telescopes of the VLTI, the field-of- 
view of the AMBER observations was only ~60 mas. Thus, 
the background contamination of the AMBER data can be ex- 
pected to be much weaker, if not negligible, compared to the 
VINCI measurements. To check this, we first performed a fit 
of the Hillier et al. (2001) model, which not only contains 
the size scaling as a free parameter, but also a fully resolved 
background component. As we expected, we found the best fit 
(smallest^ 2 ) with no background contamination. 2 Therefore, 
when we finally compared the AMBER observations with the 
model from Hillier et al. (2001), we did not introduce a back- 
ground component. In Fig. 6 (second row, left) we plot both 
the AMBER visibilities (no background correction required) 
plus the background-corrected VINCI data (assuming a 55% 
background contribution; blue triangles). As can be seen from 
the figure, these VINCI points nicely match the AMBER data 
and the corresponding fit of the NLTE model from Hillier 
et al. (2001). Therefore, from the analysis of the continuum 



2 An additional argument in favor of only a very faint background 
contribution in the AMBER UT observation can be found in the shape 
of the high spectral resolution line: the light from the speckle objects 
B, C, and D is produced in areas with velocities smaller than 50 km/s. 
Therefore, it produces a narrow emission line which should appear in 
the center of the broad Bry line. Just looking at the shape of the line, 
it can be concluded that such an effect is negligible. 
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data, we can conclude that the background contamination in 
the AMBER measurements is negligible and that the AMBER 
measurements are in good agreement with both the previous 
VINCI measurements and the model predictions from Hillier 
et al. (2001). 

3.3. Elongated shape of the continuum intensity 
distribution 

To look for detectable elongations of the continuum intensity 
distribution, we fitted an elliptically stretched 2-D version of 
the radiative transfer model visibilities from Hillier et al. (2001) 
to the measured visibilities. Our best^ 2 fit reveals a projected 
axis ratio of £ = 1.18 + 0.10 and PA = 120 ± 15°. Comparison 
with the results found by van Boekel et al. (2003) shows that the 
projected axis ratio ^ derived from the AMBER data is in basic 
agreement with the ^-broad-band values of £ = 1.25+0.05 
and PA = 138 + 7° from van Boekel et al. (2003). 

We also studied the elongation inside the Bry emission 
line at A — 2. 166 /mi, following the same procedure as in the 
continuum; i.e., we fitted an elliptically stretched 2-D version 
of the Hillier et al. model shown in Fig. 6 to the AMBER 
data. However, since the global shape of the model function 
at A = 2. 166 /mi shows stronger deviations from the measure- 
ments than in the continuum, the elongation determination suf- 
fers from larger uncertainties, resulting in large error bars of 
the fit parameters. For instance, for A = 2.166 /mi we obtained 
f = 1.66 ± 0.60 and PA=81+40° from the best ellipse fit. 

The 2-D ellipse fitting was also performed for the contin- 
uum near the He I emission line and in the center of the He I 
line (A = 2.057 /mi), where our model fits give an axis ratio 
of £ = 1.35 ± 0.30 and a PA of the major axis of 98 ± 40° in 
the continuum, and £, = 1 .74 + 0.60 and PA = 159 + 40° in the 
center of the He I emission line. It should be noted that for the 
He I line region, only four visibility points are available, cov- 
ering the small PA range of only 60°. Because of this limited 
number of data points and the small PA coverage, we conclude 
that the He I elongation measurements in the continuum as well 
as the line region are not reliable and abandoned in the further 
elongation analysis of the He I data. 

^.Frorn the #-band VINCI data, van Boekel et al. (2003) de- 
rived a PA of 138+7° for the major axis, very well aligned with 
the Homunculus (132°, Davidson et al. 2001) and in agreement 
with our results (PA = 120 ± 15°). Van Boekel's and our con- 
tinuum elongation measurements favor the physical model ac- 
cording to which r\ Car exhibits an enhanced mass loss in po- 
lar direction as proposed, for instance, by Owocki et al. (1996, 
1998) or Maeder & Desjacques (2001) for stars rotating close 
to their critical rotation speed. Axis ratios of the order of 1 .2 
appear reasonable in the context of such polar-wind models. 
Suppose, for example, that the wind's polar/equatorial density 
ratio is 2 at any given radius r, as reported by Smith et al. (2003) 
to explain latitude-dependent changes in the Balmer line pro- 
files. Relevant absorption and scattering coefficients have ra- 
dial dependencies between n e ~ r~ 2 (Thomson scattering) and 
n 2 ~ r~ A (most forms of thermal absorption and emission). 
A meridional map of projected optical thickness through the 



wind would show cross-sections of prolate spheroids, corre- 
lated with the appearance of the configuration. With the ra- 
dial dependencies and polar/equatorial density ratio mentioned 
above, these spheroids have axial ratios between about 1 .2 and 
1.4; i.e., appreciably less than 2. Viewed from an inclination 
angle i * 45°(Davidson et al. 2001), the apparent (projected) 
axis ratios are between 1.1 and 1.2. This is merely one exam- 
ple, and we have omitted many details, but it illustrates that the 
polar/equatorial density ratio is around 2, in agreement with 
Smith et al. (2003). 

Finally, Smith et al. (2003) suggested that the stellar wind 
should become basically spherical during an event at perias- 
tron. This prediction can be tested if VLTI/ AMBER data are 
obtained at the next periastron passage. 

3.4. Continuum-corrected visibilities 

3.4.1. Continuum-corrected visibility in the 
Bry emission line 

To investigate the brightness distribution in the Bry line in 
more detail, we tried to disentangle the continuum and pure line 
emission from both the AMBER data as well as the model data 
to derive the size of the pure Bry line-emitting region. Since 
the visibility measured inside an emission or absorption line is 
the composite of a pure line component and an underlying con- 
tinuum, the measured line visibility, V to t (see Fig. 7 top), has to 
be corrected for the continuum contribution to obtain the visi- 
bility Viine of the line emitting (absorbing) region. As discussed 
in Malbet et al. (2006), VW can be calculated if the continuum 
level within the line is known. If the continuum within the line 
is equal to the continuum level outside the line for an optically 
thin environment, as assumed in Fig. 7, we obtain: 

r T ^tot ' ^*tot — ^cont " ^cont / , s 

Hine = UJ 

r line 

with F tot = F cont + Fn ne being the total measured flux and 
V to t being the measured visibility (also see Fig. 7 top for illus- 
tration). As outlined in the Appendix, taking a non-zero differ- 
ential phase <t>' into account introduces an additional term in 
Eq. (1), leading to 

coat | 2 - 2V tot F tot Vcont^cont ' COS <£' 

Hine = ( 

r line 

We applied Eq. (2) to the line-dominated AMBER data 
in the 2. 155-2. 175 /mi Bry wavelength range to derive the 
continuum-corrected visibility of the region emitting the Bry 
line radiation. One big uncertainty in this correction is the un- 
known continuum flux within the line. Due to intrinsic absorp- 
tion, the continuum flux might be considerably lower within the 
line than measured at wavelengths outside the line. Especially 
the blue-shifted wing of the Bry emission line might be af- 
fected by the P Cygni-like absorption, as discussed below. In 
case such an absorption component is present, our continuum- 
correction would overestimate the size of the line-emitting re- 
gion. The influence of this effect on the red-shifted wing of the 
line is likely to be much smaller, as P Cygni-like absorption 
mainly affects the blue-shifted emission. 
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Fig. 7. Top: The solid and dashed red lines show the continuum- 
corrected visibility inside the Bry line of the HR measurement with 
the shortest projected baseline (29 m). The continuum correction fol- 
lows Eq. (2). Except for the phase, all quantities entering Eq. (2) are 
indicated in the figure. We assume F tot = F cont + F\ mc and that F cont 
inside the line is approximately equal to the level outside the line. The 
continuum-corrected visibility in the blue-shifted region of the emis- 
sion line is shown with a dashed line since it is highly uncertain due 
to the presence of the P Cygni-like absorption component (see text 
for details). Bottom: Continuum-corrected AMBER visibilities (filled 
blue squares) in the red region of the Bry line as a function of spatial 
frequency. To derive the visibilities according to Eq. (2), the data in 
the wavelength range 2.1661-2.1670 yum (red-shifted line region) were 
averaged before the continuum correction. The solid red line shows 
the continuum-corrected model visibility. Obviously, observations at 
shorter baselines are needed to further test the model predictions. 



The visibility across the Bry line is shown in Fig. 7 (top) 
for the HR data corresponding to the shortest projected base- 
line. The results for the other data sets are similar. ^From Fig. 7 
(top), one can see that after the subtraction of a continuum con- 
tribution equal to the continuum outside the line, the visibility 
reaches very small values in the center of the emission line. 
This means that the pure line-emitting region is much larger 
than the region providing the continuum flux. 

Fig. 7 (top) shows a strong asymmetry between the blue- 
and red-shifted part of the visibility in the line with respect to 
the spectrum. While the visibility (V tot as well as Vn ne ) rises 



concomitantly with the drop of the line flux on the red side, the 
situation is very different on the blue side line center. In agree- 
ment with the model predictions from Hillier et al. (2001), this 
indicates the existence of a P Cygni-like absorption component 
in this wavelength region. In fact, at A = 2.1625 jum, we see 
a small dip in the Bry spectrum in both the model spectrum 
from Hillier et al. (2001) and the HR AMBER observations. 
If such an absorption component is present, it can explain the 
asymmetric behaviour of the line visibility with respect to the 
spectrum. The P Cygni absorption in the blue wing of the Bry 
line makes the continuum correction of the visibility uncertain 
for wavelengths shorter than the central wavelength A c of the 
emission line. Because of this uncertainty, in Fig. 7 (top) the 
continuum corrected visibility is shown with a dashed line for 
A < A c , and the following discussion is restricted to the red- 
shifted region of the Bry line emission. 

The continuum-corrected AMBER visibilities in the red- 
shifted region of the Bry line are displayed in Fig. 7 (bottom) 
for all data sets. To derive the visibilities in the red region, 
the data in the wavelength range 2. 1661-2. 1670/^m were av- 
eraged before the continuum correction. To now compare the 
continuum-corrected AMBER visibilities with the model pre- 
dictions (2. 166 1-2. 1670 ^m), we constructed a model intensity 
profile of the pure Bry emission line region by subtracting the 
Hillier et al. intensity profile of the nearby continuum from the 
combined line + continuum profile. 

As Fig. 7 illustrates, the model prediction is in agreement 
with the low visibilities found for spatial frequencies beyond 
60 cycles/arcsec. On the other hand, the figure also clearly in- 
dicates that measurements at smaller projected baselines are 
needed to further constrain the Hillier et al. model in the line- 
emitting region. With the baseline coverage provided by the 
current AMBER measurements, we obtain a FWHM diameter 
of > 15.4 mas (lower limit) for the (continuum-corrected) line- 
emitting region in the red line wing. 

3.4.2. Continuum-corrected visibility in the He I 
emission line 

As can be seen in Fig. 4, the AMBER spectrum of the He I line 
shows a P Cygni-like profile with a prominent absorption and 
emission component. This is in agreement with earlier findings 
by Smith (2002) from long-slit spectroscopy using OSIRIS on 
the CTIO 4m telescope. 

To estimate the spatial scale of the region emitting the He I 
emission line, we followed the same approach as outlined in the 
previous section for the Bry line; i.e., we first applied Eq. (2) 
and then compared the continuum-corrected visibility with the 
continuum-corrected radiative transfer model of Hillier et al. 
(2001). Figure 8 (top) shows the measured flux and visibility 
for the MR-2004- 12-26 measurement with the shortest pro- 
jected baseline (43m) as well as the continuum-corrected vis- 
ibility across the He I emission component (solid red line). 
Because of the P Cygni-like absorption component, the con- 
tinuum subtraction is highly uncertain in the blue region of the 
emission line (dashed red line in Fig. 8), as already discussed 
in the context of the Bry line in Sect. 3.4.1. 
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In Fig. 8 (bottom), the continuum-corrected visibility of all 
AMBER data in the red region of the He I emission line (aver- 
aged over the wavelength range 2.057-2.058 /mi) is shown as 
a function of spatial frequency. As the figure reveals, similar to 
the Bry emission, the visibilities inside the He I emission line 
region reach rather low values. As the comparison shows, the 
line visibilities predicted by the model are much higher than the 
line visibilities measured with AMBER, indicating that the size 
of the line-emitting region in the model is too small. Rescaling 
of the model size by a factor of 2.4 results in a much better 
agreement between the model and observations (green curve in 
Fig. 8, bottom) and a FWHM diameter of > 8.2 mas, which 
is 3.6 times larger than the FWHM diameter of 2.3 mas in the 
continuum. Due to the lack of interferometric data at small pro- 
jected baselines, this value can only give a rough lower limit of 
the size. 

The results for the visibility inside the He I line can possi- 
bly be explained in a qualitative way in the framework of the 
binary model for the central object in 77 Car (e.g. Davidson et al. 
1999; Davidson 2001; Pittard & Corcoran 2002; Hillier et al. 
2006; Nielsen et al. 2006). In a model of this type, He I emis- 
sion should arise near the wind-wind interaction zone between 
the binary components. The hot secondary star is expected to 
ionize helium in a zone in the dense primary wind, adjoining 
the wind-wind interaction region. Such a region can produce 
He I recombination emission. 3 The wind-wind shocked gas, by 
contrast, is too hot for this purpose, while the density of the fast 
secondary wind is too low. Since the AMBER measurements 
(Dec. 2004 and Feb. 2005, at orbital phases <f> = 0.268 and 
<p = 0.299, see Table 1) were obtained at an intermediate phase 
between periastron in July 2003 and apoastron in April 2006, 
the extension of the He I emission zone is expected to be rather 
diffuse and larger than the continuum size. In other words, the 
He I emission zone should be fairly extended and larger than 
the Hillier et al. model prediction, which is in agreement with 
the AMBER data. 

3.5. Differential Phases and Closure Phases 

The measurement of phase information is essential for the re- 
construction of images from interferometric data, but such an 
image reconstruction is only possible with an appropriate cov- 
erage of the uv plane. Nevertheless, even single phase measure- 
ments, in particular of the closure phase and differential phase, 
provide important information. 

The closure phase (CP) is an excellent measure for asym- 
metries in the object brightness distribution. In our AMBER 
measurements, as illustrated in Figs. 2 and 4, we find that the 
CP in the continuum is zero within the errors for all the vari- 
ous projected baselines of the UT2-UT3-UT4 baseline triplet, 
indicating a point-symmetric continuum object. However, in 

3 For a qualitative sketch of the geometry, see "zone 4" in Figure 8 
of Martin et al. (2006), even though this figure was drawn to repre- 
sent He ++ in a different context. For reasonable densities, the predicted 
He + zone has a quasi-paraboloidal morphology. In addition, some ex- 
tremely dense cooled gas, labeled "zone 6" in the same figure, may 
also produce He I emission. 
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Fig. 8. Visibility in the He I line. Top: The figure is similar to Fig. 7 
(top) but displays the MR-2004- 12-26 data with the shortest projected 
baseline (43m). As in the case of the Bry emission, the figure reveals 
that the region dominated by line emission is fully resolved by the 
AMBER measurements. As in Fig. 7 (top), the continuum-corrected 
visibility is shown with a dashed line in the blue-shifted region of the 
emission line to indicate that in this region, the line visibility is highly 
uncertain due to the presence of the P Cygni absorption. Bottom: 
Continuum-corrected AMBER visibilities (filled blue squares) in the 
red region of the He I emission line as a function of spatial frequency 
and continuum-corrected visibility curve (solid red line) according to 
the NLTE model of Hillier et al. (2001). To derive the visibilities ac- 
cording to Eq. (2), the data in the wavelength range 2.057-2.058 /jm 
were averaged before the continuum correction. The model curve pre- 
dicts much larger visibility values in the line, which indicates that the 
size of the line-emitting region in the model is probably underesti- 
mated by the model. A rescaling of the size of the model by a factor 
of ~2.3 (green line) reveals better agreement with the AMBER mea- 
surements. 



the line emission, we detect a non-zero CP signal in all data 
sets. In both MR measurements covering the Bry line, we find 
the strongest CP signal in the blue wing of the emission line 
at A = 2.164 /mi (-34° and -20°) and a slightly weaker CP 
signal in the red wing of the emission line at A = 2. 167 /mi 
(+12° and +18°). We also detected non-zero CP signals in the 
HR measurement around Bry taken at a different epoch. In the 
case of the He I line, a non-zero CP could only be detected at 
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A — 2.055 pm, just in the middle between the emission and 
absorption part of the P Cygni line profile. 

The differential phase (DP) at a certain wavelength bin 
is measured relative to the phase at all wavelength bins. 
Therefore, the DP measured within a wavelength bin contain- 
ing line emission yields approximately the Fourier phase of the 
combined object (continuum plus line emission) measured rel- 
ative to the continuum. This Fourier phase might contain con- 
tributions from both the object phase of the combined object 
and a shift phase, which corresponds to the shift of the photo- 
center of the combined object relative to the photocenter of the 
continuum object. Significant non-zero DPs were detected in 
the Doppler-broadened line wings of the Bry line. Particularly 
within the blue-shifted wings, we found a strong signal (up to 
60°), whereas the signals are much weaker within the red- 
shifted line wings. These DPs might correspond to small pho- 
tocenter shifts, possibly arising if the outer Bry wind region 
consists of many clumps which are distributed asymmetrically. 
The small differential phases of up to ~ -15° for the differ- 
ent baselines of the blue-shifted light in the He I line can per- 
haps also be explained by the above-mentioned asymmetries or 
within the framework of the binary model discussed in previ- 
ous sections. In the binary model, a large fraction of the He I is 
possibly emitted from the wind-wind collision zone, which is 
located between the primary and the secondary (Davidson et al. 
1999; Davidson 2001; Pittard & Corcoran 2002; Hillier et al. 
2006; Nielsen et al. 2006). 

3.6. Modeling with an inclined aspherical wind 
geometry 

The goal of the modeling presented in this section is to find a 
model which is able to explain several remarkable features in 
our data; in particular, (a) the asymmetry in the Bry line profile 
(showing less emission in the blue-shifted wing than in the red- 
shifted wing) and the P Cygni-like absorption dip in the blue- 
shifted Bry wing, (b) the strong DP in the blue-shifted wing 
and a weaker DP signal in the red-shifted wing, and (c) the 
structure of the CP, showing a change in the sign between the 
blue- and red-shifted line wing. We aimed for a geometrical 
but physically motivated model which would reproduce these 
features at all wavelength channels simultaneously. For this, we 
concentrate on the Bry line, as this line shows a stronger phase 
signal than the He I line and was measured with a better uv 
coverage. 

As Smith et al. (2003) convincingly showed, the stellar 
wind from r\ Car seems to be strongly latitude-dependent, 
with the highest mass flux and velocities at the poles. This 
anisotropy can be understood in the context of theoretical mod- 
els (see, e.g., Maeder & Desjacques 2001), which take the 
higher temperatures at the poles (g e ff-effect) and the equatorial 
gravity darkening on a rapidly rotating star into account (von 
Zeipel effect, Zeipel 1925). As these models are quite suc- 
cessful in explaining the bipolar structure of the Homunculus 
nebula, we investigated whether such bipolar geometries with 
a latitude-dependent velocity distribution might also be suited 
to explain our interferometric data. 



Due to its success in reproducing both the spectrum and the 
measured visibilities, we based our wind model on the spheri- 
cal Hillier et al. (2001) model and superposed a weak aspher- 
ical stellar wind geometry, which is inclined with respect to 
the line-of-sight. Our model includes three components (see 
Fig. 9); namely, 

(1) a continuum component (using the Hillier et al. continuum 
CLV, see Fig. 6 top) with a blue-shifted absorption compo- 
nent, 

(2) a spherical stellar wind (using the Hillier et al. continuum- 
subtracted Bry CLV), and 

(3) an aspherical wind geometry, represented by a 4 1 ° inclined 
ellipsoid. 

The relative contribution of these different constituents to the 
total flux is given by the input spectra shown in the upper row 
of Fig. 10. For the spherical and aspherical wind component, 
we assume Gaussian-shaped spectra. The original Hillier et al. 
CLVs slightly underestimate the size of the observed structures 
(see Fig. 5). Therefore, we rescaled them by 10% to obtain a 
better agreement for the visibilities at continuum wavelengths. 

The aspherical wind of r\ Car is simulated as an ellip- 
soid with an inclination similar to the inclination angle of the 
Homunculus (41°, Smith 2006). While the south-eastern pole 
(which is inclined towards the observer) is in sight, the north- 
western pole is obscured. The latitude-dependent velocity dis- 
tribution expected for the rj Car wind was included in our model 
by coupling the latitude-dependent brightness distribution of 
the ellipsoid to the wind velocity. At the highest blue-shifted 
velocities, mainly the south-eastern polar region contributes to 
the emission (see Fig. 9a). In the red-shifted line wing, mainly 
the (obscured) north-western pole radiates (see Fig. 9b). The 
axis of the ellipsoid was assumed to be oriented along the 
Homunculus polar axis (PA 132°, Smith 2006) and its axis ratio 
was fixed to 1.5. 

As our simulations show, such an asymmetric geometry can 
already explain the measured DPs and CPs with a rather small 
contribution of the asymmetric structure to the total flux (see 
black line in Fig. 10, upper row). Although the large number 
of free parameters prevented us from scanning the whole pa- 
rameter space, we found reasonable agreement with a size of 
the ellipsoid major axis of 8 mas. Fig. 10 shows the spectrum, 
visibilities, DPs, and CPs computed from the model. 

As our model was inspired by physical models, but does not 
take the complicated radiation transport and hydrodynamics in- 
volved in reality into account, we would like to note that our 
model allows us to check for consistency between the consid- 
ered geometry and the AMBER spectro-interferometric data, 
but can neither constrain the precise parameters of a possible 
aspherical latitude-dependent stellar wind around 77 Car, nor 
can it rule out other geometries. We summarize some quali- 
tative properties of our wind model as follows: 

- The strong underlying spherical component mainly ac- 
counts for the very low visibilities measured within the line. 

- The aspherical wind component introduces the asymmetry 
required to roughly explain the measured phase signals. In 
particular, it reproduces the larger DPs and CPs within the 
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a) Blue-shifted line wing 

A = 2.1635 yum 



b) Red-shifted line wing 

A = 2. 1669 fim 



Br;' spherical wind Continuum 

emitting region 

Br; latitude-dependent, 
asphericaJ stellar wind 




Fig. 9. Left: Illustration of the 
components of our geometric 
model for an optically thick, 
latitude-dependent wind (see text 
for details). For the weak aspher- 
ical wind component, we draw 
the lines of latitudes to illustrate 
the 3D-orientation of the ellip- 
soid. Right (a, b): The upper 
row shows the brightness distri- 
bution of the modeled aspheri- 
cal wind component (item (3) in 
the text) for two representative 
wavelengths. The figures below 
show the total brightness distri- 
bution after adding the contribu- 
tions from the two spherical con- 
situents of our model. 




blue-shifted line wing compared to the red-shifted wing, as 
the red-shifted emission region is considerably obscured. It 
also accounts for the flip in the CP sign, as the photocenter 
of the line emission shifts its location between the blue- and 
red-shifted wing relative to the continuum photocenter. 
- The absorption component which we introduced in the 
blue-shifted wing of the Bry line allows us to reproduce 
the asymmetry measured in the shape of the Bry emission 
line profile (showing an increase of flux towards red-shifted 
wavelengths) and the weak dip observed at far-blue-shifted 
wavelengths. Furthermore, with the decrease of the contin- 
uum contribution, the absorption component helps to lower 
the visibilities in the blue-shifted line wing, simultaneously 
increasing the asymmetry in the brightness distribution (in- 
creasing the phase signals). Finally, with the interplay be- 
tween the absorption and emission component, our simula- 
tion reproduced a "bump" in the visibility similar to the one 
observed on the shortest baseline of our HR measurement 
(A * 2.1631ml). 



3.7. Feasibility of the detection of the hypothetical hot 
companion and the wind-wind interaction zone 

One of the most intriguing questions regarding rj Car is whether 
or not its central object is a binary, as suggested to explain cycle 
(e.g., Damineli 1996). 

3.7.1 . A simple binary continuum model 

To investigate whether the AMBER measurements presented 
here can shed more light on the binarity hypothesis, we used 
the following approach: We constructed a simple binary model 
consisting of a primary wind component with a CLV according 
to the continuum model of Hillier et al. (2006; FWHM diam- 
eter ~ 2.3 mas; see upper panels in Fig. 6) and an unresolved 
binary companion represented by a point-like source (uniform 
disk with ~ 0.1 mas FWHM diameter). The secondary com- 
ponent is predicted to be approximately located at PA -36° 
with a separation of 8 mas from the primary for the time of 
the AMBER observations (Nielsen et al. 2006). The continuum 
flux ratio q was treated as a free parameter. We would like to 
note that in our model, we assumed that all /T-band light from 
the secondary is reaching us unprocessed; i.e. we ignored a pos- 
sible dilution or re-distribution of the secondary's radiation. 
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Fig. 10. Observables computed from our optically thick, latitude-dependent wind model (see Fig. 9 for a model illustration). The points 
(crosses) represent the measurements (as also shown in Fig. 2), and the solid lines give the observables computed from our model. The upper 
row shows the contributions from the various model components to the total flux. Besides the continuum emission (purple line), we introduced 
a spherical (blue line) and an aspherical (black line) wind component. 
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a) Continuum emission with binary 




c) 2D- Visibility with binary 
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Fig. 11. Simulation illustrating the signatures of a binary companion at the predicted position (for the orbital phase at the time of our continuum 

observations around the Bry line; separation ~ 8 mas, PA 36° ; see Nielsen et al. 2006). For the primary component, we assume a continuum 

emission CLV from Hillier et al. (2006), to which we add a uniform disk of 0. 1 mas diameter to account for the binary companion (see a). 
Figure b shows the 2D-visibility function of a single component model containing only the continuum-emitting object with Hillier-type CLV. 
Figure c shows the 2D- visibility corresponding to the binary model shown in panel a for a small A"-band flux ratio q = 10 between the 
continuum emitting region (Hillier-type CLV) and the companion. In d the residuals of the visibilities (upper panel) and closure phases (lower 
panel) between the models with and without a companion are shown as a function of flux ratio q for baselines and PAs corresponding to the 
seven AMBER measurements. The dashed horizontal lines indicate the uncertainty of the measurements. Given these uncertainties, the figure 
illustrates that a companion signature should be visibile from the AMBER closure phase measurements, if the binary is less than ~ 110 times 
fainter compared to the primary. See text for details. 



We calculated the 2D visibility function (see Fig. 11c) of 
this model intensity distribution for different values of q, as 
well as the closure phases for the baselines and PAs corre- 
sponding to our AMBER measurements (Fig. 1 Id). Finally, we 
compared the results with those obtained from a single compo- 
nent model where only the primary wind is present (Fig. 1 lb). 
The differences of the visibilities and closure phases between 
the single star and the binary model (at the baselines and PAs 
corresponding to our AMBER measurements) are displayed in 



Fig. lid as a function of the K-band flux ratio of the binary 
components. 

Fig. lid shows two interesting results: First, the closure 
phase is more sensitive to the binary signature than the vis- 
ibilities and, thus, a more suitable observable to constrain the 
binary hypothesis. And second, given the accuracies of our first 
AMBER visibility and closure phase measurements (indicated 
by the horizontal dashed-dotted lines), we can conclude, for 
the particular model shown in Fig. 11, that the AMBER clo- 
sure phases put an intensity ratio limit q^ n a 1 10 on the binary 
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TT-band flux ratio. This limit is in line with the estimate q « 200 
given by Hillier et al. (2006). Thus, based on the model shown 
in Fig. 11, the AMBER measurements are not in conflict with 
recent model predictions for the binary. 

To investigate whether we can put similar constraints on the 
minimum TT-band flux ratio q m \ n for arbitrary separations and 
PAs, we calculated a larger grid of binary models and compared 
the residuals of visibilities and closure phases analogue to the 
example shown in Fig. 1 1 . 

For these grid calculations, we used values in the range 
from 4 to 14 mas for the binary separation with increments of 
1 mas, and PAs of the secondary in the whole range from 0° 
to 360° in steps of 10°. The /iT-band flux ratio q of the binary 
components was varied in the range from 1 to 250 with Aq = 4. 
As a result of the grid calculation, we obtained the minimum 
/f-band flux ratio g m ; n as a function of binary separation and 
orientation. 

Whereas the study with a fixed companion position pre- 
sented above allowed us to put rather stringent constraints on 
<7min (see Fig. 11), this systematic study revealed that due to 
the rather poor uv coverage, a few very specific binary pa- 
rameter sets exist where we are only sensitive to q m \ n « 10. 
Nevertheless, for the above-mentioned separation interval (4 to 
14 mas), we found that we are able to detect companions up 
to q m in = 50 at more than 90% of all PAs. In order to push 
this sensitivity limit in future observations, a better uv cover- 
age will be required. Together with the expected higher closure 
phase accuracy, AMBER will be sensitive up to q > 200 and, 
therefore, have the potential to probe the currently favored bi- 
nary models. 

3.7.2. Can AMBER detect a He I wind-wind interaction 
zone shifted a few mas from the primary wind? 

In the context of the binary hypothesis, it is also important to 
discuss the implications for the interpretation of the AMBER 
He I measurements. According to the binary model, a large 
fraction of the He I line emission should arise from the wind- 
wind collision zone expected between the primary and the sec- 
ondary (Davidson et al. 1999; Davidson 2001; Hillier et al. 
2006; Nielsen et al. 2006). The exact intensity ratio of primary 
He I wind and He I emission from the wind-wind interaction 
zone is not known. Figure 5 suggests that during the AMBER 
observations, the total He I flux was roughly two times larger 
than the model prediction of Hillier et al. (2001) for the primary 
He I wind. 

At the orbital phases of the AMBER measurements, the 
wind-wind collision zone should be at resolvable distances 
from r] Car's primary (resolution ~ 5 mas; companion sepa- 
ration ~ 8 mas, PA ~ -36°, Nielsen et al. 2006). Looking at the 
AMBER He I data, we see that the differential as well as the 
closure phases are zero everywhere except for the transition re- 
gion between the absorption and emission part of the He I line, 
where we find differential phases of ~ 10-20° and a closure 
phase of ~-30°; i.e., the phases measured across the He I line 
are significantly weaker compared to the Bry line. The question 
is now, why AMBER measured weaker phase signals within 



the He I line and if this result is in line with the predictions of 
the wind-wind collision model. 

One possible explanation for the small measured phases 
could be the orientation of the binary orbit. If the orbit's ma- 
jor axis is nearly aligned with the line-of-sight, the photocen- 
ter shift inside the He I line will be very small. In addition, 
the deviations from point symmetry would be rather small. 
Therefore, in the case of this special geometry, both differen- 
tial phases and the closure phase would be small, in qualitative 
agreement with the AMBER data. Another explanation could 
be that the contribution of the wind-wind collision zone to the 
He I line emission is much weaker than that of the primary 
wind. However, this is not very likely (see Hillier et al. 2006). 

A different explanation for the weak phases can be found 
from a modeling approach similar to the one for the Bry line 
region outlined in Sect. 3.6. Based on the results presented in 
Sects. 3.1 and 3.4.2, we constructed a simple He I model con- 
sisting of a spherical primary wind component with a Hillier- 
type CLV (2.5 mas FWHM diameter) and an extended spher- 
ical He I line-emitting region with Gaussian CLV and a 7 mas 
FWHM diameter (i.e., for simplicity, we assumed that all He I 
flux is emitted from the wind-wind interaction region; however, 
some fraction of He I is also emitted from the primary wind; see 
Hillier et al. 2006 and Fig. 5 of the present paper). The center 
of the line-emitting component of this model is located 3 mas 
away from the primary wind component towards PA 132°; i.e., 
in the direction of the Homunculus axis. The spectra of the con- 
tinuum and line-emitting components were chosen in such a 
way that the combined spectrum resembles the observed He I 
line spectrum. 

The modeling results show that this simple model is ap- 
proximately able to simultaneously reproduce the observed 
spectrum and the wavelength dependence of visibilities, dif- 
ferential phases (10-20°), and closure phases ( — 30°). Thus, 
our simple model example illustrates that the AMBER mea- 
surements can be understood in the context of a binary model 
for r\ Car and the predicted He I wind-wind collision scenario 
(e.g. Davidson et al. 1999; Davidson 2001; Hillier et al. 2006; 
Nielsen et al. 2006). We note that the model parameter val- 
ues given above are of preliminary nature. A more detailed, 
quantitative modeling is in preparation and will be subject of a 
forthcoming paper. Furthermore, we would like to emphasize, 
as already discussed in previous sections, that there are likely 
to be three sources of He I emission - the primary wind, a wind- 
wind interaction zone (bow shock), and the ionized wind zone 
caused by the ionization of the secondary. For both the bow 
shock and the ionized wind zone, the ionizing UV radiation 
field of the secondary is of crucial importance. On the basis of 
the observed blue-shift and the weakness of the He I during the 
event, we believe that the primary wind contribution is small. 
It is not yet possible to decide on the relative contributions of 
the bow shock and the ionized wind region. 

4. Conclusions 

In this paper we present the first near-infrared spectro- 
interferometry of the enigmatic Luminous Blue Variable 77 Car 
obtained with AMBER, the 3-telescope beam combiner of 
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ESO's VLTI. In total, three measurements with spectral res- 
olutions of R = 1,500 and R = 12,000 were carried out in 
Dec. 2004 (cf> = 0.268) and Feb. 2005 ((f> = 0.299), covering 
two spectral windows around the He I and Bry emission lines 
at A = 2.059 and 2.166/im, respectively. From the measure- 
ments, we obtained spectra, visibilities, differential visibilities, 
differential phases, and closure phases. From the analysis of the 
data, we derived the following conclusions: 

- In the /iT-band continuum, we resolved 77 Car's optically 
thick wind. From a Gaussian fit of the /iT-band continuum 
visibilities in the projected baseline range from 28-89 m, 
we obtained a FWHM diameter of 4.0 + 0.2 mas. Taking 
the different fields-of-view into account, we found good 
agreement between the AMBER measurements and previ- 
ous VLTI/VINCI observations of 77 Car presented by van 
Boekel et al. (2003). 

- When comparing the AMBER continuum visibilities with 
the NLTE radiative transfer model from Hillier et al. 
(2001), we find very good agreement between the model 
and observations. The best fit was obtained with a slightly 
rescaled version of the original Hillier et al. model (rescal- 
ing by 1-2%), corresponding to FWHM diameters of 2.27 
mas at A — 2.040 /um and 2.33 mas at A = 2. 174 /mi. 

- If we fit Hillier et al. (2001) model visibilities to the ob- 
served AMBER visibilities, we obtain, for example, 50% 
encircled-energy diameters of 4.2, 6.5, and 9.6 mas in the 
2.17/um continuum, the He I, and the Bry emission lines, 
respectively. 

- In the continuum around the Bry line, we found an asym- 
metry towards position angle PA=120° + 15° with a pro- 
jected axis ratio of £ = 1.18 + 0.10. This result con- 
firms the earlier finding of van Boekel et al. (2003) using 
VLTI/VINCI and supports theoretical studies which predict 
an enhanced mass loss in polar direction for massive stars 
rotating close to their critical rotation rate (e.g. Owocki 
et al. 1996, 1998). 

- For both the Bry and the He I emission lines, we measured 
non-zero differential phases and non-zero closure phases 
within the emission lines, indicating a complex, asymmet- 
ric object structure. 

- We presented a physically motivated model which shows 
that the asymmetries (DPs and CPs) measured within the 
wings of the Bry line are consistent with the geome- 
try expected for an aspherical, latitude-dependent stellar 
wind. Additional VLTI/ AMBER measurements and radia- 
tive transfer modeling will be required to determine the pre- 
cise parameters of such an inclined aspherical wind. 

- Using a simple binary model, we finally looked for a pos- 
sible binary signature in the AMBER closures phases. For 
separations in the range from 4 to 14 mas and arbitrary PAs, 
our simple model reveals a minimum TT-band flux ratio of 
-50 with a 90% likelihood. 

Our observations demonstrate the potential of 
VLTI/AMBER observations to unveil new structures of 
77 Car on the scales of milliarcseconds. Repeated observations 
will allow us to trace changes in observed morphology over 
77 Car's spectroscopic 5.5 yr period, possibly revealing the 



motion of the wind-wind collision zone as predicted by the 
77 Car binary model. Furthermore, future AMBER observations 
with higher accuracy might be sensitive enough to directly 
detect the hypothetical hot companion. 
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Appendix A: Wavelength calibration 

To obtain both an accurate wavelength calibration of the 
AMBER raw data and properly calibrated spectra of rj Car, we 
compared the AMBER raw spectra of r\ Car as well as the cal- 
ibrator stars L Car, HD 93030, and HD 89682 with a /ST-band 
telluric spectrum recorded at the Kitt Peak Observatory with 
a spectral resolution of 40,000. For the comparison with the 
AMBER spectrum, this telluric spectrum was spectrally con- 
volved to match the spectral resolution of the AMBER mea- 
surements with high (R ~ 12,000) and medium (R ~ 1,500) 
spectral resolution. 

The result of the comparison is shown in Fig. A. 1 . In the 
upper panel, the high spectral resolution AMBER spectra of 
j] Car and the calibrator L Car are shown together with the tel- 
luric spectrum with R = 10, 000. ^From the comparison with 
the telluric spectrum, we identified 7 prominent telluric absorp- 
tion features in the L Car spectrum, which are indicated by the 
dashed vertical lines. The strongest absorption line seen in the 
L Car spectrum is not telluric, but can be identified as intrinsic 
Bry absorption in L Car. Therefore, to properly calibrate the 
i] Car spectrum with the L Car spectrum, we had to interpolate 
the Bry line region in the L Car spectrum before dividing the 
two spectra. From the spectral calibration shown in Fig. A.l, 
we estimated a wavelength calibration error of the AMBER 
data AA = 3 x 10" V m 

The lower panel in Fig. A.l shows the wavelength calibra- 
tion of the medium spectral resolution data in the wavelength 
region around the He I line. The figure contains the two rj Car 
MR spectra and the spectra of the two corresponding calibrator 
stars, HD 93030 and HD 89682, as well as the telluric spectra 
with spectral resolutions of R = 40, 000 and R = 1, 500. As the 
telluric spectra reveal, there is a forest of telluric lines in the 
spectral region around the He I line. As can be seen in Fig. 4, 
the modulation of the continuum flux introduced by the telluric 
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Fig. A.l. Top: Spectral calibration of the AMBER observations. The 
figure shows the spectrum of r\ Car (blue) and the calibrator star 
L Car (green) around the Bry line obtained from the HR measure- 
ments in Feb. 2005. In addition, a telluric spectrum from the Kitt Peak 
Observatory is shown in red. The original telluric spectrum with a 
spectral resolution of R = 40, 000 was spectrally convolved to match 
the resolution of the AMBER measurements. For the flux calibration 
of the rj Car spectrum, the Bry line in the L Car spectrum was inter- 
polated before the division. See text for further details. Bottom: Same 
as top panel, but for both MR measurements around the He I emis- 
sion line. Together with the rj Car spectra, the spectra of the calibra- 
tors HD 93030 and HD 89682 are displayed. In addition, the telluric 
spectrum obtained at Kitt Peak with R = 40, 000 and a spectrally con- 
volved telluric spectrum is shown, which matches the spectral resolu- 
tion of the MR measurements (R = 1, 500). The figure illustrates that 
the forest of telluric lines forms a quasi-continuum which modulates 
the AMBER spectra. 



quasi-continuum cancels out completely when the rj Car spec- 
tra are divided by the corresponding calibrator spectra, which 
show no prominent intrinsic line features. Since there are no 
sharp spectral features in the 2.03-2.08 /im region of either the 
calibrator or telluric spectras which could be used for the spec- 
tral calibration, we estimated a wavelength calibration error 
AA = 6 x 10~ 4 yum for the MR He I data. On the other hand, for 
the MR data around the Bry line, we found AA - 4-x 10~ 4 yum. 



Appendix B: Continuum uniform disk and Gauss 
diameter fits 

For each spectral channel as well as for an averaged contin- 
uum, we performed 1-D fits to the visibility data using simple 
uniform disk (UD) and Gaussian models. In this step of the 
analysis, possible asymmetries were ignored and all visibility 
points at a given wavelength were fitted together, regardless of 
the position angle of the observations. The results of these 1-D 
fits are illustrated in the two upper panels of Fig. B. 1 for the av- 
eraged continuum data in the wavelength ranges 2.03-2.08 /im 
and 2.155-2.175 /im, respectively. As the figure reveals, neither 
a uniform disk nor a single Gaussian provides a good fit to the 
continuum data. At least, this is true as long as no contamina- 
tion by a fully resolved background component is taken into 
account. 
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Fig. B.2. Dependence of the Gaussian FWHM diameter on the fit 
range. The figure shows the background-corrected visibilities obtained 
with VLTI/VINCI (see van Boekel et al. 2003) as well as Gaussian fits 
of (a) all four data points (long-dashed green line), (b) only the point 
with q = 45 cycles/arcsec (short-dashed blue), (c) only the point cor- 
responding to the longest baseline (dotted purple), and (d) only the 
point corresponding to the shortest baseline (dashed-dotted light blue). 
See the labels for the Gaussian FWHM diameters resulting from the 
different fits. The figure illustrates that the fitted diameter strongly de- 
pends on the spatial frequency range which is used to fit the data. The 
strong diameter variation (in this case, the diameter changes by a fac- 
tor of ~ 3) occurs since a Gaussian is not a good representation of the 
measured visibility function. See text for further discussion. 



The wavelength dependence of the apparent size obtained 
from the UD and GAUSS model fits for the individual spectral 
channels is shown in the lower left panel of Fig. B. 1 . This panel 
illustrates that the equivalent UD and GAUSS /iT-band diame- 
ters of rj Car derived from the AMBER data are ~4 and ~6.5 
mas, respectively. 

It should be added here that a good fit of the AMBER 
data using, for instance, a Gaussian can indeed be obtained 
when a certain amount of contamination due to a fully re- 
solved background component is taken into account (see also 
Petrov et al. 2006a). To illustrate that, we performed Gaussian 
fits to the AMBER data, where we introduced such a fully 
resolved component as a free fitting parameter. We found 
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Fig. B.l. 1-D visiblity fits of the AMBER continuum mea- 
surements. The upper two panels show uniform-disk and 
Gaussian fits for the averaged continuum around the Hel 
and Bry lines, respectively. As both panels illustrate, neither 
a uniform disk nor a Gaussian is a good representation of the 
AMBER measurements. The panel on the left summarizes 
the fitting attempt for all continuum wavelength channels 
by displaying the fitted diameters of the UD and GAUSS 
models as a function of wavelength. 



that the best Gaussian fit is obtained with a FWHM diameter 
^Gauss = 3.01 mas and a ~ 29% background contamination 
for the Hel continuum region and c/Gauss = 3.32 mas and a 
~ 30% background contamination for the Bry continuum re- 
gion. Thus, from this fit we would derive a background con- 
tamination which is only ~ 50% smaller than in the VINCI 
data. We think that such a large amount of background con- 
tamination is not very likely given the small AMBER fiber 
aperture (60 mas) of the 8.2 m telescopes. We think that the 
large amount of background contamination needed to find a 
reasonable Gauss fit just reflects the fact that a Gaussian is not 
appropriate to describe the observations. This is confirmed by 
the fact that for the fit of the radiative transfer model of Hillier 
et al. (2001), no background component has to be taken into 
account to reproduce the AMBER measurements. 

We would like to note here that the Gaussian FWHM diam- 
eters of typically ~ 4 mas found from the AMBER measure- 
ments are not in contrast to the value dc- dU s>i ~ 7 mas found by 
van Boekel et al. (2003) from VLTI/VINCI observations for the 
following reason: Since a Gaussian is not a good representation 
of both the VINCI and the AMBER visibilities, the diameter 
resulting from a Gaussian fit strongly depends on the fit range. 
This is illustrated in Fig. B.2 for the four VINCI measurements 
given in Fig. 1 of van Boekel et al. (2003). As the figure shows, 
from a Gaussian fit of all four data points, c/Gauss = 6.0 mas is 
obtained. If only the data point with q = 45 cycles/arcsec is 
fitted (corresponding to a projected baseline length of ~ 24 m), 
we get o?Gauss = 7.65 mas. This is in agreement with the val- 
ues given in van Boekel et al. (2003) for the elliptical Gaussian 
fit of the large number of VINCI measurements with a pro- 
jected baseline of 24 m (see their Fig. 2). On the other hand, if 



we fit only the VINCI data point corresponding to the longest 
projected baseline (q = 136 cycles/arcsec), a Gaussian fit pro- 
vides a?Gauss = 4.39 mas (see Fig. B.2), which is very close to 
the diameter we obtain from the AMBER measurements for 
A ~ 2.174/im ((^Gauss = 4.35 mas). This is not surprising since 
the spatial frequency of this VINCI data point agrees with the 
average spatial frequency of our AMBER observations {q ~50- 
200 cycles/arcsec). Thus, it can be concluded that good agree- 
ment between the Gaussian FWHM diameters derived from the 
AMBER and VINCI measurements is found if a comparable 
spatial frequency range is used for the fit. 



Appendix C: Visibility and differential phase of an 
emission line object 

We assume that the target's intensity distribution can be 
described by two components: the continuum spectrum 
o C ont(x,y, X) and the emission line spectrum on ne (x,y; A). In the 
part of the spectrum containing the emission line, both o cont and 
oii n e contribute to the total intensity distribution o tot . According 
to the van-Zittert-Zernike theorem, the Fourier transforms 
O cont (B/A) and On ne (B/A) of o cont (x,y; A) and o hne (x,y;A) are 
measured with an optical long baseline interferometer at wave- 
length A and projected baseline vector B. In the following, we 
assume that all Fourier spectra are normalized to 1 at frequency 
zero. The complex Fourier spectrum O tot (B/A) of the intensity 
distribution o tot measured at the emission line An ne is given by 

Otot = — (^cont ■ (5 CO nt + ^line - <5line), (C.l) 

r cont ~t" r line 
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where F cont and Fi; ne are the fluxes of the continuum com- 
ponent o cont and the line component oii ne , respectively. In the 
emission line, the total flux measured is F tot := F cont + Fu ne . 

^From the spectrally dispersed interferometric data, we can 
derive the differential phase, which is the difference of the 
Fourier phases of the continuum component o cont and the to- 
tal intensity o tot in the emission line. The differential phase 
<J>'(B//tii ne ) in the emission line at A\i ne is given by 

Ocont ' 0* tot = V cont ■ V t ot ' , (C2) 

where Q>'(B/A Me ) := Ow^Miine) - ® tot (B/ A line ). 
^contiB / An ne ) is the Fourier phase of the continuum com- 
ponent, and Oto^B/zliine) denotes the Fourier phase of o tot at 
the position of the emission line An ne . The asterisk * in this 
equation denotes conjugate complex operation. V mnt (B / An ne ) 
describes the visibility of the continuum component at the 
position of the emission line Ai\ ne , and Vtot(fiMiine) is the 
visibility measured at the position of the emission line Anne- 
Inserting Eq. (C.l) into Eq. (C.2) yields 

Vrnnt ■ Vtot ■ e 



(Fcont Vcont + Fine Vline<?' A *), (C.3) 



+ Fat 



where A<&(B I A\\ ne ) denotes the difference of the Fourier phases 
of the continuum and line components; i.e., AO(B / A\ me ) := 

3>cont(5Mline) ~ ^line^Mine)- <J>cont(#Mline) and OlineCBMine) 

are the Fourier phases of the continuum and line components, 
respectively. 

In the vector representation of complex numbers, the three 
quantities F C ontVcont, fWVrme, and F tot V tat form a triangle 
with one corner placed at the center of the coordinate system. 
According to the law of cosines, the correlated flux of the line 
component is given by Eq. (2) (see Sect. 3.4.1): 

IFine Vlinel 2 = l^totVtot| 2 + l^cont V C ont| 2 ~ 

~2 ■ f^ot^ot • fcont Vcont ■ cos(<D'). (C.4) 

Since the flux F\i ne can be calculated from the measured fluxes 
Fcont and F tot , the visibility Vij ne of the line component can be 
derived using Eq. (C.4). 

Applying the law of sines to this triangle in the complex 
plane yields the differential phase A<$>(B / A\ me ), which is the dif- 
ference between the Fourier phase O CO nt(B/^iine) of the contin- 
uum component and the Fourier phase <t>ii ne (B//lii ne ) of the line 
component: 

sin(A(D) = sin(<D')--^lM, (C.5) 

Inline Hine I 

where <&'(B/A]i„ e ) is the differential phase measured at the po- 
sition of the emission line An De . 



